Limits on Cosmological Variation of Strong Interaction and Quark Masses 
from Big Bang Nucleosynthesis, Cosmic, Laboratory and Oklo Data 
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Recent data on cosmological variation of the electromagnetic fine structure constant from distant 
quasar (QSO) absorption spectra have inspired a more general discussion of possible variation of 
other constants. We discuss variation of strong scale and quark masses. We derive the limits on 
their relative change from (i) primordial Big-Bang Nucleosynthesis (BBN); (ii) Oklo natural nuclear 
reactor, (iii) quasar absorption spectra, and (iv) laboratory measurements of hyperfine intervals. 
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I. INTRODUCTION 

Time variation of major constants of physics is an old 
and fascinating topic, its discussion by many great physi- 
cists - Dirac as the most famous example - had surfaced 
many times in the past. Recent attention to this issue 
was caused by astronomical data which seem to suggest 
a variation of electromagnetic a at the 10~^ level for the 
time scale 10 bn years, see The issue discussed in 
this work is related to it, although indirectly. Instead of 
looking into atomic spectra and testing a stability of elec- 
tric charge, we will discuss possible variations of nuclear 
properties induced by a change in strong and weak scales. 
We will not go into theoretical discussion of why such 
changes may occur and how they can be related to mod- 
ification of electromagnetic a. Our aim is to identify the 
most stringent phenomenological limitations on such a 
change, at (i) a time of the order of few minutes, when 
the Big Bang Nucleosynthesis (BBN) took place, as well 
as (ii) at the time of Oklo natural nuclear reactor (1.8 bn 
years ago) , (iii) when quasar radiation has been absorbed 
in the most distant gas clouds (3-10 bn years ago) and 
(iv) at the present time. 

Mentioning relevant literature we start with the BBN 
limits on electromagnetic a, obtained in . The main re- 
sults come from variation of late-time nuclear reactions. 
Because of low temperatures and velocities involved at 
this stage, those reactions have quite significant suppres- 
sion due to Coulomb barriers, in spite of the fact that 
only Z=l-3 is involved. These limits arc in the following 
range 
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In general, all the models for time variations of electro- 
magnetic/ weak/strong interactions can be divided into 
two distinct classes, depending on whether it originates 
in (i) infrared or (ii) ultraviolet. The former approach 
ascribe variations to some hypothetical interaction of the 
corresponding gauge bosons with some matter in Uni- 
verse, such as vacuum expectation values (VEVs) or 
"condensates" of some scalar fields. Those typically have 
zero momentum but can have cosmological time depen- 
dence. We would not discuss it: for recent example and 
references see [||. 



We would however mention few details from two re- 
cent examples of the latter approach, by Calmet and 
H. Fritzsch ||] and Langacker, Segre and Strassler 0. 
Their main assumption is that Grand Unification j^] of 
electromagnetic, weak and strong forces holds at any 
time. Therefore, a relation between all three coupling 
constants exists: truly modified 2 parameters are in this 
approach the unification scal^AouT and the value of the 
unified coupling acuT at this scale. Their time variation 
is assumed to propagate down the scales by the usual 
(unmodified) renormalization group. 

If this assumption is correct, any variation of elec- 
tromagnetic a should be accompanied by a variation of 
strong and weak couplings as well. Specific predictions 
need a model, we will mention the one discussed in [Q. 
In it the QCD scale Aqcd (determined as usual by a 
continuation of the running coupling constant into its - 
unphysical - Landau pole) is modified as follows 
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Another focus of our work is possible limits on cos- 
mological modifications of quark masses. According to 
Standard Model, they are related to electroweak symme- 
try breaking scale, as well as to some Yukawa couplings 
hi . In 1^ running of those has been considered, with 
a (model-dependent) conclusion that quark mass indeed 
may have a different (and stronger) change 
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Large coefficients in these expressions are generic for 
GUT and other approaches, in which modifications come 
from high scales: they appear because weak and strong 
couplings run more. 



* One might think that if the GUT scale be used to set units, 
its variation would be impossible to detect without explicit 
measurements related to gravity. But it is not so, since the 
cosmological expansion itself (which is quite important for 
BBN) contains the Newton's constant (or the Plank mass) in 
the Hubble constant. 
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If a coefHcients of such a magnitude are indeed there, 
at the BBN time of few minutes the QCD scale and quark 
masses would be modified quite a bit, if the upper limit 
(|l|) be used in the r.h.s. 

The type of questions we are trying to answer in this 
work are: Do we know whether it might or might not ac- 
tually happened? Which simultaneous change of strong 
and weak interaction scales is or is not observable? What 
observables are the most useful ones, for that purpose? 
What are the actual limits on their variation which can 
be determined from BBN and other cosmological and lab- 
oratory data? 

Let us repeat, that although we use the above men- 
tioned papers as a motivation, we do not rely on any 
particular model. Nevertheless, we will at the end of the 
paper return to these predictions in order to see whether 
our limits on time variation imply stronger or weaker ef- 
fects than the electromagnetic ones. 

II. THE ROLE OF HEAVY AND STRANGE 
QUARK MASSES IN HADRONIC/NUCLEAR 
OBSERVABLES 

Both papers just mentioned argue for what we 
would call the zeroth approximation to QCD modifica- 
tion. It assumes that the QCD scale Aqcd is so dominant 
in all hadronic and nuclear phenomena that all dimen- 
sional parameters - hadronic masses, magnetic moments, 
energies of nuclear levels, etc - are to a good approxima- 
tion simply proportional to its respective powers. 

If so, any time variation of the overall strong inter- 
action scale would not change dimensionless quantities 
(such as mass ratios or g-factors) which we can observe. 
The absolute scale of hadronic and nuclear spectroscopy 
may change, but even if we would be able to observe it 
from a cosmological distances, such a modification would 
easily be confused with the overall redshift. 

Fortunately, this pessimistic situation]^ is in fact rather 
far from reality. Quark masses do play significant role in 
hadronic/ nuclear physics, and, if they have time mod- 
ification different from that of Aqcd, those can be de- 
tected. 

Logistically it is convenient to start with masses of 
heavy - c,b,t - quarks. They do play a role in running of 
the strong charge, from high scales down, changing beta 
function each time a corresponding mass scale is passed. 
However, as it is well known, those effects can be readily 
absorbed in re-definition of Aqcd- 

If strong coupling is the same at some high normaliza- 
tion point M, between nib and nit, the corresponding re- 



' Motivated historically by large number of colors limit or 
quenched lattice QCD. Both may be reasonable starting ap- 
proximations, which are not expected to be really accurate. 



lations between Aqcd with all experimental quarks and 
without heavy c,b,t quarks is as follows: 

Awithout~c,b,t = Aqcd{ ■ )^^^'' (4) 

Note that rather small powers of the masses are in- 
volved in this relation. In effect, we indeed may pre- 
tend that c,b,t quarks do not exist at all, as far as basic 
hadronic/nuclear physics is concerned. 

The situation is completely different with the next 
quark flavor we have to discuss, the strange quark. It 
is still true, that if one fixes strong coupling as(fc) at 
some sufficiently high scal^ and then considers the role 
of non-zero in perturbative beta function, the effect 
is negligible. (Since its scale is too low for QCD to be 
passed by.) 

But hadronic/nuclear masses and properties are not 
determined by perturbative diagrams, leading to beta 
function: they are of course determined by much more 
complicated non — perturbative dynamics. Although it 
is far from being completely understood, it is clear that 
it does indeed depend strongly on quark masses. In par- 
ticular, the "strange part" of the vacuum energy density 
can be estimated, because the derivative of the vacuum 
energy 

r>f 

^ =< 0|ss|0 -1.4 /m-3 (5) 

is known||. Thus the linear term in the strange part of 
the vacuum energy 

e,=ms< ss >« -0.2 GeV/fm^ (6) 

is not negligible compared to gluonic vacuum energy |^ 

^ -(0.5-1) GeV/fm^ (8) 

making 20-40 percent of it. (The numerator in this ex- 
pression is the familiar coefficient of the QCD beta func- 
tion, with Nc = 3, Nf = 3 being the number of colors 
and relevant flavors. It appears because this expression, 
known as the scale anomaly to lowest order, has the same 
origin as the beta function itself.) 



'It is in fact done on the lattice, where k is the inverse lattice 
spacing, typically 2-3 GeV. 

^For definiteness, this number had come from the QCD sum 
rules, which typically correspond to operator normalization at 
fj, = 1 GeV. The lattice numbers are similar, but normalized 
at inverse lattice spacing, typically fi = 2 GeV. Anomalous 
dimension of this operator lead to small difference between 
two normalizations which we ignore here. 
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Furthermore, for the nucleon one finds that similar 
"strange fractions" of their masses are of the same mag- 
nitude, e.g.Q 



dm 



N 



drUs 



=< N\ss\N >« 1.5 



(9) 



Putting into hnear expansion the strange quark mass 
rus = 120 - 140 MeV one finds that about 1/5 of the 
nucleon mass comes from the "strange term" . 

(Although in this paper we cannot go into discussion of 
why it is the case, let us make a small digression. First of 
all, the reader should not be confused with the fact that 
only a very small fraction of the energy of a fast moving 
nucleon is due to the strange sea, as experimentally mea- 
sured partonic densities tell us. Such drastic difference 
between vector-like (or chiral-even) and scalar-like (or 
chiral-odd) operators is very common feature of the non- 
perturbative QCD. Its origin is related with dominance 
of the instanton-induced t Hooft interaction, see |jl^. In 
short, it happens because topological tunneling events 
should necessarily involve all light fermion flavors.) 

Returning to cosmology, we conclude that due to 
strange terms, any variation of quark masses would im- 
ply significant modification of hadronic masses and other 
properties. It remains a challenging task for model- 
builders and lattice practitioners to establish to what 
extent the O(rns) part of hadronic masses is or is not 
universal. 

In principle, what we would call the most pessimistic 
scenario is possible, in which Aqcd and enter into 
all hadronic observables in one combination 



Ae// = A 
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where K is some universal constant. If such a scenario 
happens to be true, its time modification can indeed be 
neutralized by a change of units, since experiments can 
only measure dimensionless ratios. 

In fact, when lattice practitioners express the obtained 
results in terms of the so called "physical units"|^, the 
dependence on indeed tends to become weaker. 

However at the moment this is just a hint, with accu- 
racy not better than say 10 percent, and there is no rea- 
son to expect this scenario to be the case. We mentioned 



such a pessimistic case provocatively, emphasizing that 
at the moment we lack solid theory which would explain 
how any particular hadronic observable depends on mg. 
In general, variation ofrug a/o77.e can noticeably influence 
strong interaction parameters, since different quantities 
in general depends differently on it. Let us give an ex- 
ample. 

One of the most important quantity for astronomi- 
cal and laboratory experiments is the magnetic moments 
of nuclei. For example, the ratio of hyperfine splitting 
to molecular rotational intervals is proportional to a^g 
where g is defined by the magnetic moment 



2mr,c 



(11) 



In zeroth approximation as well as in the most pessimistic 
scenario discussed above, only one dimensional parame- 
ter {Aqcd and Aqcd + Knig, respectively) exist, so a 
dimensionless g factor cannot have any time variation. 
But, we repeat, there is no general argument for such ap- 
proximations to be accurate. The magnetic moment and 
the nucleon mass are not directly related to each other 
at the QCD IcvclF^, and the dimensionless derivative 



(12) 



remains unknown even for the proton and neutron, to say 
nothing about the composite nuclei. If it is not the same 
as derivative (^ for the nucleon mass, any time variation 
of JTis/A would induce a variation of the nuclear g factors. 

The role of light quark masses is another issue, and 
a part of magnetic moments related to the contribution 
of the so called "pion cloud" will be briefly discussed in 
section 



XI 



In summary, any dimensionless ratios should be viewed 
as a function of the ratio 
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which for small variations is reduced to partial derivatives 
such as mentioned above. 



III. THE ROLE OF LIGHT QUARK MASSES 



** The reader may find discussion of phenomenological sit- 
uation in tiie first paper from [0, while the second contains 
lattice calculations of this quantity. For reference, their con- 
clusion the r.h.s. of (^) is 1.53± 0.07, with the errors being 
statistical only. 

t+Those units are defined by some non-perturbative observ- 
able such as (i) p- meson mass or (ii) the string tension, or 
(iii) the force between two point-like charges at fixed distance. 
All of those are measured on the lattice, with whatever quark 
masses one wants to have, and then put equal to its observable 
value in real world by a decree. 



Unlike strange quark mass, we have more solid the- 
ory explaining what the effect of a change of light quark 



**For example, in non-relativistic quark model the nucleon 
mass is approximately 3 times constituent quark mass, and 
quark magnetic moment is given by "quark magneton" . How- 
ever even in this model there is also a binding energy and other 
corrections. A constituent quark itself is a complicated com- 
posite object, so one should not expect its magnetic moment 
to be exactly equal to the Dirac value related to its mass. 
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masses niu, rud relative to that of the strong scale can be. 
The main focus of this work is this particular variation, 
which at the end we will be able to constrain rather well. 

Since the pion is a Goldstone boson, its mass scales as 
a geometric mean between weak and strong scales]^ 

2 



(14) 



Therefore the appropriate parameter characterizing the 
relative change of the pion mass ratio to the strong scale 
can be defined as 



A, 
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(15) 



i5(^^)/(^!^) 

2 ^QCD ^QCD 

Because the pion mass determines the range of nuclear 
forces, its modification leads directly to changes in nu- 
clear properties. The main question we would like to 
study is how such a change in the pion mass relative 
to that of other hadrons, described by a non-zero S^r, is 
limited at the cosmological time when primordial nucle- 
osynthesis took place. 

The first limit on relative change of quark masses have 
been put in M, 



-0.1 < 6{- 



< 0.02 



(16) 



where both numerical values come from the observational 
uncertainty of He* production, and T^, is the freezeout 
temperature for neutrino-induced weak processes. This 
effect is however not very restrictive, for the following 
reasons: 

(i) The sensitivity of He* production to any variation is 
itself not very impressive. (Below we will discuss d or 
Li^ yields which may vary by orders of magnitude and is 
much more sensitive.); 

(ii) As correctly explained in [Q, T^, scales as 



'/Ml 



(17) 



as a function of weak scale (Higgs VEV) v and Plank 
mass, while to„ — rup variation is expected to come pri- 
marily from S^rud — ?ti„), with smaller electromagnetic 
correction. Therefore, both numerator and denominator 
in (^6|) mostly reflect the same physics and thus a large 
portion of the modifications, if exist, would tend to cancel 
in this combination. 



The QCD anomalous dimension of the quark mass is can- 
celled by the opposite one of the quark condensate, in the 
GOR relation: thus quark masses here are meant with the 
dependence on the normalization point removed. The same 
remark should also be made about the term with the strange 
quark. (We thank T. Dent, who reminded us that this com- 
ment is needed.) 



IV. CRUDE LIMITS ON SCALE VARIATIONS AT 
THE BBN TIME 

Before we come to specifics, let us explain our ba- 
sic philosophy in selection of the observables. To get 
the maximal sensitivity, one has to focus on phenomena 
which may vary by very large factors. Basically, there are 
two sources for that, inside the dynamics which drives the 
BBN. One is the Boltzmann factor^ which may at late 
times reach 10 orders of magnitude or more; the other 
is Gamow factor due to Coulomb barriers, which may 
reach 3-4 orders of magnitude. So naturally, production 
of the heaviest primordial nuclides - especially Li^, which 
is sensitive to both of those - would be most promising 
place to look. 

We start however with a preliminary discussion of pos- 
sible drastic changes when (5^ is of the order of several 
percents, and then return to more delicate limits on scale 
variation. 

The nucleosynthesis starts with the two-nucleon states. 
Standard nuclear forces produce one bound state - the 
deuteron, pn with T=0,J=1 - and a virtual triplet of 
states - pp,pn,nn with T=1,S=L=J=0. 

If pion was lighter at the BBN time relative to its 
present value, it leads to better binding. Rather dra- 
matic effects should have occured if the virtual states 
unbound by standard nuclear forces were bound. 

Specifically, a binding of (pp) state, with its subse- 
quent beta decay into the deuteron, could eliminate free 
protons, in obvious contradiction with observations [^ . 
Conditions for a binding of pn and nn states have been 
studied recently in |22|. These states may add new paths 
to nucleosynthesis. 

Note also, that if the pion was so much heavier than 
now that the deuteron gets unbound, no primordial nu- 
cleosynthesis could possible proceed at all. We will not 
discuss those in this paper, as they are superseded by the 
modification of the deuteron discussed below. 

As the only exception to this rule, we briefly discuss 
another important bottleneck on the way toward heavy 
nuclides, the absence of the A=:5 bound states. We have 
shown below that with (J^r ~< —0.052 it can be bridged 
in a modified world. 



V. DEFINING TWO MOST IMPORTANT 
TEMPERATURES OF THE BBN 

As we will argue throughout this work that the most 
sensitive parameter in BBN remains the binding energy 
of the deuteron \Ed\ {Ed is negative), we provide a (very 
brief) account of BBN, with emphasis on the role of 
deuterons. For details of the evolution, the role of all 
processes involved the reader should consult appropriate 
reviews, e.g. 

When the basic reaction producing deuterons p + n ^ 
d -f 7 is in equilibrium at high enough T, the density of 
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deuterons relative to photons is of the order of 



nd 



(18) 



Here 77 « 3 * 10"^" is the famous primordial baryon-to- 
photon ratio. Although it remains unknown what have 
created it and even at which stage of the Bing Bang it 
happened, we will assume it happens early and do not 
consider its variations. 

So, the density of deuterium remains negligible small 
till the Boltzmann factor (needed for the photons to split 
d) helps. In particular, when this factor is so large than it 
can compensate one of the 77, the deuteron fraction may 
be comparable to that of p,n. This condition determine 
the first crucial temperature value, to be denoted by Td 



Tjexp{\Ed\/Td) ^ 1 



(19) 



Note that Td will be modified below together with the 
value of Ed- At standard BBN parameters it is about 70 
keV, at which Ud reaches its maximum, see fig. 1(a). 

What happens after d reaches this maximum, is sig- 
nificant reduction of the density of d (and other species) 
due to several reactions leading into the best bound light 
nuclei, He'^. At such low T the equilibrium configuration 
would wipe out all of other nuclides. Indeed, an advantage 
in binding of iJe* by more than 20 MeV the Boltzmann 
factor is enormous. 

However, the universe expansion does not allow to 
reach such equilibrium. Its rate leads to the final freeze- 
out of all production reactions, the stage when reaction 
rates and the Hubble expansion rates are comparable. 
The Hubble rate, according to one of the Friedman eqn. 
for fiat Universe and zero cosmological constant (not im- 
portant so early anyway) is 



= [R/Rf ^ (87r/3)G'Are 



(20) 



where Gjq , e are Newton's constant and mater energy 
density. Ignoring numerical constants, H ~ T'^/Mp, 
with Mp being the Plank mass, and ignoring the T- 
dependence of the reaction rate itself < av >t, one can 
obtain the freeze-out temperature 



Tf ~ l/{r]Mp < av >Tf) 



(21) 



Note how small parameter rj fight with large Mp. For 
standard BBN it is Tf « 35 keV. 



VI. HOW VARIATION OF THE DEUTERON 
BINDING AFFECT THE BBN 



identically, so that this change can be eliminated from 
the discussion by simply using modified units. In such 
units all reactions rates, bindings and Td would be the 
same as in the standard BBN. The only difference ap- 
pears in the freeze-out time and temperature T/. The 
reason is, as we changed the units to adjust to time vary- 
ing scales, the value of the Plank mass entering relation 
(pi| ) gets modified instead. 

We can estimate crudely the effect of that modification 
as follows. At the second stage of the BBN, Td <T <Tf, 
a decrease of d and Li^ is roughly power-like 



fd 



Li7 



T- 



a = 6-7 



(22) 



The fall is by about 2 orders of magnitude. Using this 
trend we conclude that the modification of their yields 
within a factor 2 (the magnitude of current error bars) 
corresponds to certain change in Tf, which implies the 
following limits on variation of the strong scale (relative 
to gravity) 



SjAQCp/Mp) 
{Aqcd/Mp 



< 0.1 



(23) 



We now switch to discussion of the relative change, be- 
tween quark masses and hadronic scale, leading to mod- 
ification of the pion-induced forces and consequently, of 
the deuteron binding energy. We will look at modifica- 
tion of Td, which is so important for BBN final observable 
yields. 
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Now we can proceed to discussion of how modification 
of the fundamental interactions would affect the BBN 
yields. 

Let us start with the most pessimistic case when 
both the strong scale and quark masses are modified 



FIG. 1. (a) Schematic dependence of the deuterium mass 
fraction fd on temperature T (MeV). (b) Schematic depen- 
dence of the deuterium and He'^ mass fractions on the binding 
energy of the deuteron Ed (MeV). The vertical dotted line in- 
dicate its experimental value, -2.2 MeV. 
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A qualitative dependence of the yields of different 
species on the magnitude of the deuteron binding is 
schematically shown in fig.l (b). As Td is allowed to 
vary, the first thing to note is that there exist an opti- 
mum for d production corresponding to the case when 
two crucial temperatures are close to each other, 



(24) 



It is easy to see from expressions given above, that it 
happens if the deuteron binding is around Ed ~ —1 MeV, 
and its fraction reaches at this time about a percent level. 

When the value of Td reduces further, so that the re- 
lation between the two temperatures is inverted and Tf 
becomes smaller than Td, the stage at which d and other 
light nuclei are "eaten" by i/e'* is no longer present. All 
heavier nuclides, t,He^, He'^ etc. reach their final yields 
from below, basically tracing the deuteron production. 
There is no time for He* to grab most of the available 
neutrons any more, and its final yield start declining (see 
Fig. 1(b)). Finally, when deuterons become nearly un- 
bound, all of the primordial nucleosynthesis is wiped out 
altogether. 

As data tell us that d fraction is at the level 10^* — 
10~^, the BBN happened clearly away from the maxi- 
mum. Furthermore, as He'^ fraction is large, we are def- 
initely at the left branch of the curve in Fig. 1(b). This 
shows exponential growth, and thus sensitivity of the re- 
sults to the variation of Ed and Td can be estimated from 
its slope, which we obtained from We conclude from 
this estimate that a change of the d yield by the fac- 
tor 2 corresponds to a constrain on relative variation of 
deuteron binding by 

(^). < 0.075 (25) 



VII. FURTHER LIMITATIONS AND THE 
COULOMB BARRIERS 

Now we can discuss the refinement of qualitative pic- 
ture described above, in which the absolute change of 
strong scale has been eliminated by a change of units, 
and the only weak effect included has been a quark mass. 
Now we let the electromagnetic effects come into the 
game, in their simplest and most important (exponen- 
tial) form. We mean the Coulomb barriers, described be 
well known Gamow factors 

2Tre^ZiZ2 , li^e^Z^Z^. 

jGamow = exp{ ) (26) 



For standard BBN we estimated that the product 
of Gamow factors of reactions leading to Li^ is about 
exp{lA * vq/v), where vq is the unmodified value***. So, 
assuming that experimental data on Li^ are restricted 
within a factor 2, we see that only a variation of 



(-^)Li7 < U.l 



(27) 



can be tolerated. It leads to a limit on S-rr comparable to 
the one obtained above from deuteron yield. 



VIII. THE MODIFIED DEUTERON 

In this section we try to relate the change of deuteron 
binding SEd to the modification of the fundamental pa- 
rameter introduced above. 

In general, it is a very non-trivial dynamical issue, 
which is far from being really understood. We have dis- 
cussed above to which extent the QCD vacuum energy 
(^) and the nucleon mass (||) depends on quark masses. 
May be one day such information will be available from 
lattice QCD for nuclear forces as well, but right now we 
do not have it and have to rely on model-dependent po- 
tentials. 

Of course, one can identify single pion exchange forces. 
Especially for the deuteron channel, those lead to well 
known tensor forces producing deuteron quadrupole mo- 
ment. A textbook 1-pion exchange corresponds to the 
following potential 



T/ f ^ exp{-m^r) 

r 



(28) 



where = 0.08. The pion mass in denominator has 
a tilde: it indicates that this mass has been put there 
by hand for normalization purposes, and, unlike masses 
in other places, it will not have any variations. (It has 
been put there in order to cancel the pion mass squared 
coming from differentiation of the exponent, which we 
will evaluate first.) 

When the variation of the pion mass is sufhciently 
small, the variation is simply given in the first order 
of the perturbation theory SEd —< 0\6V\0 > For any 
weakly bound state the wave function outside the poten- 
tial range can be approximated by simple expression 



exp{—Kr)/r 



(29) 



where parameter n is related to binding energy by Ed = 
/nriM- If we use this expression till the core size, we 
can easily evaluate it 



Here v is the relative velocity, which scales as (T jm)^!'^ . 
As explained above, the relevant T varies in between Td 
and Tf. As the former one is modified, together with Ed, 
the Gamow factors change accordingly. 



***We have also checked that this estimate agrees with the 
results by Q of the dependence of its yield on a, which come 
from running a compete BBN code. 
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21^7^^; — ^ ^ — \ 

We therefore get from these terms 



5\Ed 
\EA 



-0.05 



(31) 



The sign tells us that reducing the pion mass we reduce 
binding: this counterintuitive result comes from the 
in pre-exponent. 

One more one-pion-exchange term we have ignored is 
the one in which two derivatives act on 1/r: this leads 
to delta function. Naively it does not contribute since at 
r=0 the wave function is vanishingly small due to repul- 
sive core. This can be cured by the account for a finite 
size of the nucleon. 

However the main problem with perturbative calcula- 
tion outlined above is that the 1-pion-exchange terms are 
actually several times smaller compared to phenomeno- 
logical potentials needed to reproduce NN scattering and 
the deuteron binding. Furthermore, the phenomenolog- 
ical potentials which fit scattering and d data actually 
ascribe most of the potential to 2 and even 3-pion ex- 
changes. As an example, we used well known work by 
Hamada and Johnston , which has the central poten- 
tial in the deuteron channel in the form of subsequent 
pion exchanges. We have only modified the pion mass in 
the exponents: 



V{r) = V^Y{x)[l + acY{x) + KY{xf 



(32) 



Y{x) 



(33) 



where x = r* m^r is distance in units of unmodified pion 
mass, and — — —1 [Q. The potential also has 
infinite repulsive core with the radius 0.4 fm. 

We have ignored tensor and spin-orbit potentials, as 
having minor effect on the deuteron binding, adjusted 
the overall coefficient to have the experimental binding 
energy -2.2 MeV. After that, we have modified the pion 
mass in (^2|), and calculated binding energy by solving 
the s-wave Schroedinger eqn. For small pion modifica- 
tion we need, we found good linear dependence with the 
following coefficient 



^■gQ^ethods is ref. Like us, these authors have shown 

examples producing opposite sign of the derivative and 
conclude that, strictly speaking, neither the magnitude 
nor even the sign can be definitely obtained at this time. 
Further extensive lattice and chiral perturbation theory 
studies are needed. 

If the derivative is small near zero, our arguments lose 
its weight. However, it would be very unnatural to get a 
small value of this derivative. The attractive and repul- 
sive parts of the potential have already conspired to get 
deuteron at the binding edge. One more fine tuning, ex- 
actly at the physical quark mass values (which does not 
have any special meaning from the QCD point of view) 
is very unlikely to happen. 



IX. BINDING OF He^ 

We now return to He"^ + n channel mentioned above, 
and evaluate which BBN limits on modification of the 
pion mass it will produce. 

In this channel there is resonance at energy 0.77 

MeV. If reduction of the pion mass can make it bound, 
it will drastically enhance production of Li. Without 
bound He^ BBN has to jump over A=5, e.g. by a reac- 
tion i/e^ -|- i, which is impaired by the Coulomb barrier 
as well as by very low concentration of t. 

We now study the sensitivity of He^ binding to mod- 
ifications of the nuclear potential. Before we discuss the 
calculations, let us make a couple of qualitative points. 
First of all, this case is more complicated compared to 
deuteron, because we now discussed very weakly bound 
states, and so no simple linear dependence on (5^ is ex- 
pected. Furthermore, let us point out an interesting 
quantum-mechanical distinction between low-lying levels 
with zero and non-zero angular momentum 1. The wave 
function at large r, outside of the interaction range, for 
zero energy level can easily be found from Schroedinger 
equation 



5\Ed\ 
\Ed\ 



(34) 



We have presented these two estimates as reasonable 
minimal and maximal hounds on this derivative. Pre- 
sumably the true value is somewhere in between 3 and 
-18. 

The issue has been discussed in literature. A gen- 
eral discussion of how all nuclear physics would change if 
quark masses, the number of flavors or even the number 
of colors be modified, can be found in . Probably the 
latest paper addressing deuteron binding with modern 



i/'(r,£; = 0) - r-(i+') 



(35) 



In the s-wave, 1=0, it is unacceptable because such wave 
function is non-normalizable solution. In other words, 
keeping exp{—Kr) in ( p9| ) is crucial, no matter how small 
K is. It is no longer so for 1=1: the tail of the wave 
function in question is now normalizable. In short, for 
I > the centrifugal barrier keeps particles inside the 
attractive potential, contrary to the I = states. This 
effect makes I > states more sensitive to the change of 
the potential. 

We have modeled the potential of the interaction be- 
tween a neutron and He^ in the following form 



Vir) 



1 + w{r/c)'^ 



' exp[{r - c)/A] -I- 1 



(36) 



where w=.445, A = .327 fm. If c=1.01 fm it is a good fit 
to experimental density distribution in He^ [fla. In the 
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potential we simply changed it to c = 1.01/m + l/m^, 
with variable pion mass. The depth of the potential Q 
has been tuned to reproduce the position of the above 
mentioned resonance, it gave Q=:-32.6 MeV. After that 
we start changing the pion mass till this level becomes 
bound, E=0. This happens at the magnitude of the pion 
mass modification 



> -0.052 



(37) 



If that would happen at the BBN time, the yields of Li 
would be dramatically enhanced, by orders of magnitude, 
contrary to observations. (This is why it is written as the 
inequality. ) 

This value can be compared to the natural small pa- 
rameter of the virtual He^ level Eres/V ^ 0.02. The 
needed shift is somewhat larger because, in spite of gen- 
eral argument given above, its wave function is spread to 
large r, see fig.Q. 

Note that this limit (|37| ) would become an order of 
magnitude weaker if we base our consideration on the di- 
rect effect of the 1-pion exchange eq.(^8|) - see discussion 
in the previous section. For He^ there is an additional 
suppression because the average value of the potential 
( p8| ) over closed Is shell is zero (the effect, however, ap- 
pears due to the exchange interaction and correlation cor- 
rections) . 

As a result, we conclude that the limit from the He^ 
binding cannot compete with that from the deuteron 
modification in its importance in BBN. 



The next bottleneck, effectively blocking synthesis of 
heavy elements, is the gap at A=8 nuclei, which too read- 
ily decay into two alpha-particles. Although we have not 
investigated this case, we do not expect it to beat the 
limits related to deuteron modification as well. 



X. LIMITATIONS FROM OKLO DATA 

Finally, let us deviate from our discussion of BBN to a 
related subject, namely similar limits following from data 
on natural nuclear reactor in Oklo active about 2 bn years 
ago. The most sensitive phenomenon (used previously 
for limits on the variations of the electromagnetic a) is 
disappearance of certain isotopes (especially Sm^'^^) pos- 
sessing a neutron resonance close to zero . Today the 
lowest resonance energy is only Eq — 0.0973 ± 0.0002 eV 
is larger compared to its width, so the neutron capture 
cross section a ^ ^/Eq. The data constrain the ratio of 
this cross section to the non-resonance one. It therefore 
implies 1^^^ that these data constrain the variation of the 
following ratio 5{Eq/Ei) where Ei ^ 1 MeV is a typical 



single-particle energy scale, which may be viewed as the 
energy of some 1-body "doorway" state. 

A generic expression for the level energy in terms of 
fundamental parameters of QCD can be written as fol- 
lows 



Ei Ai * Aqcd + Bi*mq + da* Aqcd 



(38) 



800 



600 



y 400- 



200- 




+A(mqAQCD)^/^ + ... 

where Ai,Bi, Ci, Di are some coefficients. The first term 
is the basic QCD term, while others are corrections due 
to quark mass, pious and electromagnetism. Without 
Bi,Di terms one can see that in the Eq/Ei ratio the 
QCD scale drops out, confirming our general statement 
above that in such kind of approximations the variation 
of Aqc d itself cannot be seen. 

The sum is very small for Eq, just because we delib- 
erately picked up the lowest resonance, but (and this is 
our main point) there is no reason to expect each term to 
be especially small. For example the electromagnetic 
term is about 1 MeV. 

Let us estimate what variation of the resonance en- 
ergy would result from a modification of the pion mass. 
As we did above for the deuteron, we assume that the 
main effect comes from increase of the radius R of the 
nuclear potential well. The energy of the resonances 
Ei = Eexoitation " Sn cousists of cxcitation energy of a 
compound nucleus, minus the neutron separation energy 
Sn- This, in turn, is a depth of the potential well V minus 



FIG. 2. The wave function for zero energy neutron in He^ 
(arbitrary normalization) versus the distance r, in fm. Note 
that it is normalizable, although it does spread to rather large 



^^^Of course, under assumption that the same resonance was 
the lowest one at the time of Oklo reactor. 
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the neutron Fermi energy ep, Sn = V — ep- The latter 
scales like 1/ if the radius of the well is changed. The 
kinetic part of the excitation energy E^xcitation scales in 
the same way. Adding both, one gets shift of the reso- 
nance 

SEi = —{ep + E excitation)— fT- (39) 

H 

For resonance near zero the combination in brackets is 
approximately V ^ MeV. If i? = 5/m + I/tott then 
6R/R = -Sm^/iRml) 

iVl =3*10'l'5.| <0.2 (40) 
t^i 

The r.h.s. above comes from the observational limits 
claimed in [ p^ . The resulting limitation on pion mod- 
ification at time « 1.8 bn years ago is 

Note that the authors of the last work in found 
also the non-zero solution = — 1 ± 0.1. This solution 
corresponds to the same resonance moved below thermal 
neutron energy. In this case 5-^ ~ — 4*10~®. In principle, 
the total number of the solutions can be very large since 
Sm}'^^ nucleus has millions of resonances and each of 
them can provide two new solutions (thermal neutron 
energy on the right tale or left tale of the resonance). 
However, these extra solutions are probably excluded by 
the measurements of the neutron capture cross-sections 
for other nuclei since no significant changes have been 
observed there also, see p^ . 

As in the cases of the deuteron and He' binding one 
may argue that the limits on 5tt presented above should 
be weaker by an order of magnitude since the direct con- 
tribution of the 1-pion exchange eq.(p8|) to the energy is 
small. To clarify this point it may be useful to perform a 
numerical calculation of this contribution to the neutron 
separation energy in Sm}^^ and He' . 

XI. LIMITS FROM ASTROPHYSICAL AND 
LABORATORY MEASUREMENTS 

Comparison of atomic H 21 cm (hyperfine) transition 
with molecular rotational transitions gave the following 
limits on y = a^gp ||l^ 

6Y/Y = (-.20 ± 0.44)10-5 for redshift z=0.2467 and 
5Y/Y = (-.16 ± 0.54)10-5 for z^O.6847. The second 
limit corresponds to roughly t=6 bn years ago. 

As we have already discussed above, only in the most 
pessimistic scenario all strong interaction phenomena de- 
pends on only one parameter, e.g. Ae// = KQCD-\-K*'rrLs, 
its time variation cannot change dimensionless quantities 
like proton magnetic g- factor gp. If so, the limits given 
above are just limits on variation of a. However, in gen- 
eral there is no reason to think this to be the case, and 



one may wander which limits can be put from these data 
on a cosmological variation of ttIs/Aqcd- 

Another issue here is a contribution proportional to 
light quark masses, mu,md. Let us first make a quali- 
tative point, suggesting that their role in magnetic mo- 
ments is expected to be lager than in hadronic masses. 
Hadrons are surrounded by the so called "pion cloud" , 
which have small virtual momenta p. Masses depends on 
it in the form \/ + m"^ while magnetic moments have 
p X f', T being 

the distance form the center. Small masses are partly 
compensated by large r in the latter but not former case. 

Model-dependent estimates support this idea. The 
magnitude of the effect varies a lot between models, 
and as an example we use rather conservative treat- 
ment by T.Sato and S.Sawada It can be seen 
as a minimal estimate: they identified the contribution 
of small virtual momenta p < A by using form-factors 
^ A^/(A^ +P^), and have shown that a consistent pic- 
ture for p,n,d and hyperon magnetic moment emerges if 
the cutoff is A ~ m^. Furthermore, this contribution is 
shown to be basically proportional to A^ ~ ~ rriq. 
They found that the cloud contribution is about 1 per- 
cent of the proton magnetic moment, but 7.3 percents for 
the neutron (to be compared to the 1 percent level for 
the masses). 

We conclude that at least due to the pion cloud effect 
one should expect that the giromagnetic ratios g for nu- 
clei have a term proportional to mught I ^qcd contribut- 
ing of the order of several percents. Combining it with 
the contribution proportional to rris/ Aqcd (which we 
hope does not exactly conspire with the effect on masses 
to cancel completely) , we expect overall effect at the level 
of 1/10 of g is of this origin. 

Assuming that it has such magnitude, that a varia- 
tion of alpha and rriq/ Kqcd do not conspire to produce 
observed zero, and a simplest linear dependence 

gp = gp{mq = Q){l + q—^) (42) 

we may interpret the above mentioned limits as the fol- 
lowing on variation of this ratio 

|5(/^)/(t^)I < 10-^ (43) 

J-^QCD i^QCD 

This should be compared with the limit on X = 
a^gprne/nip @ SX/X = (0.7 ± 1.1) * 10"^ for z=1.8. 
This limit was interpreted as a limit on variation of a or 
me/mp. It also can be viewed as a limit on variation of 
mq/ Kqcd- Although few times weaker (~ 2* 10""*) than 
the limit (^3|), it corresponds to higher redshift. 

The limits on variation of TUq/ Kqcd can also be ob- 
tained from laboratory measurements of ratios of hyper- 
fine splittings. By comparison the rates of two clocks 
based on different atoms, say H and Hg'^ ||l^, we com- 
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pare ^ g-factors of quite different nuclei . In terms of 
the standard shell model description, this gives the ratio 
of proton and neutron spin g-factors. Other examples, 
such as using hyperfine transition in Cs as a frequency 
standard, would also involve the orbital g-factor, with 
gi = 1. In principle, corrections to shell model include 
"exchange currents" which also contribute to magnetic 
moments. All of the above may have different depen- 
dence on mq/ Aqcd, and we conclude that 



S{Aqcd/Mp 
{Aqcd/Mp 



< 0.1 



(47) 



d Ai d niq 
— In — — = il — In - — — 

dt Ai dt Aqcd 



(44) 



where Ai, A2 are hyperfine structure constants of differ- 
ent atoms, where K, a combination of derivatives (^ 12), 
remains unknown but can be as big as 1/10 or even larger. 
Using such tentative value and H, Cs and Hg~^ mea- 
surements p^ , pO[ |, we obtain the limit on variation of 
mq/AqcD about 5 * 10~^^ per year. 



XII. SUMMARY 

Combining our strongest limits on the deuteron bind- 
ing, from deuteron ( [25| ) and Li^ (p7|), corresponding to 
variation of their production by a factor 2, with a relation 
between modification of the deuteron binding and modi- 
fication of the pion mass ( |3^ ) . Both effects suggest about 
the same BBN limit on the modification of the pion mass 
relative to st rong interaction scale Aqcd- 

Using eq. ( p4]) we obtain: 



\5Abbn < 0.005 
Eq. (|3l| ) provides a more conservative limit 
ISttIbbn < 0.03 



(45) 



(46) 



and we think the true limit is somewhere in between. 

We have investigated other effects, such as binding of 
He^ or pp, nn, up (S=0) states, but found that in these 
cases the needed pion modification about an order of 
magnitude larger. (It is expected, since all these states 
are more loosely bound than the deuteron.) If mg mod- 
ification relative to strong scale are as large as our limit 
on light quark modification just mentioned, it means the 
nucleon mass can be modified within ±2 MeV due to 
strange term. Note also, that our limit on quark mass 
modification is stronger than the limit [Q| coming from 
proton- neutron mass difference ([l^). 

We also pointed out significantly weaker limits on a 
simultaneous modification of strong scale and ruq scale 
at the same rate, relative to the gravity scale 



Limits on a variation of the quark mass relative to strong 
scale at the 10"'' level 3-10 bn years ago follows from 
observations of distant objects ([43|), while at the time 
1.8 bn years ago the Oklo data lead to even better limits, 
at the 10"^ - 10"^ level. 

Although there is no general relation between varia- 
tion of weak, strong and electromagnetic constants, as 
we mention in the Introduction it is implied by Grand 
Unification . If one uses those (||,^) , one finds that 
all our limits on relative weak/strong modification are 
much more restrictive than the corresponding limits on 
the modification of the electromagnetic a. In the case 
of astronomical observations, in which variation of alpha 
seems to be seen, one may either soon find the variation 
of g factors, or rule out relations between couplings based 
on Grand Unification idea. 

Finally, let us emphasize that our discussion is semi- 
qualitative in many aspects, and a lot of quantitative 
work remains to be done. Theory-wise, the most straight- 
forward thing to do is to add modifications directly into 
the BBN code, and get more quantitative limits. Al- 
though there seem to be no particular problem with the 
standard BBN at the moment, it is still true that the 
calculated yields and observations typically differ by 1-2 
standard deviations Q|. Therefore it seem to be worth 
wile to make a global fit to data with unrestricted mod- 
ification parameters (like our S,^) and see whether zero 
value would or would not be the best one. 

Another challenge to the theory, probably mostly lat- 
tice simulations, is to clarify the issue of the dependence 
of various hadronic parameters on the strange quark mass 
nis, especially how universal are the derivatives like (|^) 
for all hadrons. 

Experimental laboratory work and astronomical obser- 
vations of distant objects can significantly enhance the 
limits available today, hopefully with a non-zero effect 
eventually observed. 
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